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Abstract. We investigate the interior hyperbolic region of axisymmetric and stationary black holes 
surrounded by a matter distribution. First, we treat the corresponding initial value problem of the hyperbolic 
Einstein equations numerically in terms of a single-domain fully pseudo-spectral scheme. Thereafter, a 
rigorous mathematical approach is given, in which soliton methods are utilized to derive an explicit relation 
between the event horizon and an inner Cauchy horizon. This horizon arises as the boundary of the future 
domain of dependence of the event horizon. Our numerical studies provide strong evidence for the validity 
of the universal relation A~ = (SvrJ)^ where A'^ and A~ are the areas of event and inner Cauchy 
horizon respectively, and J denotes the angular momentum. With our analytical considerations we are able 
to prove this relation rigorously. 



1. Introduction 

Axisymmetric and stationary black hole space-times are characterized by the existence of two Killing 
vectors ^ and rj. Outside the black hole, these vectors can always be linearly combined to form a time- 
like vector. In contrast, any such non-trivial linear combination, taken in some interior neighborhood 
of the event horizon H'^ , inevitably leads to a space-like vector. The corresponding Einstein equations, 
expressed in Boyer-Lindquist-type coordinates, are elliptic in the black hole's exterior and hyperbolic in 
its interior. It has been shown in [JJ that for non- vanishing angular momentum J of the black hole, there 
always exists a regular boundary of the future domain of dependence of the event horizon, the 'inner 
Cauchy horizon' 'H~ . 

In this contribution, we investigate the interior region of axisymmetric and stationary black holes 
which are surrounded by a matter distribution. The Einstein field equations are considered in Boyer- 
Lindquist-type, as well as Weyl coordinates, where they possess a well-defined behavior at the horizons. 
Note that in Boyer-Lindquist-type coordinates, the mathematical form of the field equations is the same 

In section [2] we apply a fully pseudo-spectral method and solve the degenerate hyperbolic Einstein 
equations as an initial value problem. We expand the field quantities with respect to both spatial and time 
coordinates in terms of Chebyshev polynomials. Starting at the past Cauchy horizon we evolve the 
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data into the black hole interior up to the future Cauchy horizon 1-L . The data, found in this manner at 
T-L^ , obey the remarkable property 

A+A- = {Snjf , (1) 

where A'^ and A~ are the areas of event and inner Cauchy horizon respectively. The numerical 
solutions possess the extreme accuracy known for spectral methods when applied to elliptic problems. 
In particular, the relation ([T) was found to be valid up to 12 decimal digits. 

In section [3] we review analytical work through which ([TJ was proved rigorously ||T]. The 
corresponding techniques are based on 'Backlund transformations', a method known from soli ton 
theory. We conclude this article with a discussion of generalized and related black hole properties in 
axisymmetry. 



2. Numerical studies 

In Boyer-Lindquist-type coordinates (i?, 9, ip, t), the metric becomes singular at ii = itrh, 

= A -f^ 2 + + nsin^ e{d(p'^ - Ludtf - -{R^ - rl)dt^ , (2) 



i?2 _ y ^ 

where itrh denote the coordinate locations of the two horizons H^. The interior hyperbohc region of 
the black hole is characterized by 

e (-r-h,rh), ee[0,TT]. (3) 

Here, the metric coefficients jl,uj,u depend on R and 9 only and are regular at i? = ztrt (see EO). 
The corresponding hyperbolic Einstein equations read as follows: 

{R^ -rl)u^RR + 2Ru^R + Ufie + UfiCoi9 = 1 - ^ s\t?9 [oj'^r + '^^^ | , (4) 



{R^-r1){uj^RR + Auj^RUfi)+uj^ee = - w e(3cot 6* + 4u e) , (5) 



{R^ - rl)tl^RR + Rfi^R + ilfie = sin^^* \^^% + ^2 ) + ^^^R 



h / 

-{R^ - rl)u^j^ - Ufi {ufi + cot 6*) , (6) 

where tt := ^ In [u/u^^ and /i := ^ In (/i/uj^) with := u{R = rY^,9 = ojB Note that fi does 
not enter in the equations dH |5]l. An explicit solution of this system is given by the Kerr metric which 
describes a single rotating black hole in vacuum. 

At i? = ±rh, the metric potentials have to obey the following boundary conditions (see IHO) 

uj = constant = uj^ , 

: <! (7) 

lb 2 Th M,_R + u^ge + cot 9 = 1 — |-ti^a;^^ sin^ 9 

with the horizon angular velocities a;^. 

^ Here we choose the north pole value of ii as normalization constant. It remains positive and finite even in the degenerate 
limit rh — > 0. 
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Expressions for the black hole angular momentum J and horizon areas have been discussed in 
El 121. They are given by 



^ I " f- - ^ I -.2, , \ 3/ 



J = ^ f d^ab = —[^J u'io^Rl^^sin-'edB, (8) 
= y V^A^I-^i sin d0 = 47rn^, 



(9) 



where the 'north pole' values of u are denoted by 

n± = u{R = ±rh,6' = 0). 

Note that, for the area formula (|9}, we used the following relations: 

At the Horizons T-L^ [R = ±rh) : \/Jiu = constant 
At the Rotation axis (sin 9 = 0) : fi = u . 



(10) 



In our numerical scheme we solve the equations dD and dS). The corresponding solutions provide us 
with the relevant information to compute J and A^. 

As a first step, we express the metric coefficients u and w in terms of auxiliary potentials U and f2: 

u{R,9) = u+{e) + {R-rh)U{R,9) 

uj{R,e) = UJ+ + {R-rh)uj+{e) + {R-ri,ff2{R,e). 

In this ansatz, we take into account the boundary condition ijj{R = rh, 6) = uj'^ ■ We prescribe freely a 
set of 'initial' data {u^ {9),ijj'^ ^u}^{6)} which is then evolved from i? = rh into the black hole interior 
up to i? = — rh via the equations ^ and ^ (see figure [T]!. Note that the second boundary condition in 
([7]) fixes the first 'time' derivative n,_R(rh, 6) in terms of u^{6) and u}^{9), and hence u,ij(rh, 6) cannot 
be prescribed freely. 

We have chosen a. fully pseudo-spectral scheme B for the numerical integration of and In 
practice, we expand the functions U and f2 with respect to Chebyshev polynomials: 

n n / u\ / o 

c 

j=0 k=0 



m)„ I R \ „ 12, 

V 

j=0 k=0 



and consider the field equations on an n x n-'spectral grid' with Gauss-Lobatto gridpoints (i?j, 9^) (see 
figure [B, 



Rj=rhCos , 6lfc = 7rsin^— — {j,k = 0, . . . ,n - 1) . (13) 

n — 1 2{n — 1) 



The pseudo-spectral collocation point method provides an approximation of the values 

Ujk = U (Rj ,ek) , Qjk = ^ (Rj , Ok) (14) 
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Figure 1. Tiie interior hyperbolic black hole region with past {R = r^) and future (R = — rh) Cauchy 
horizons T-L^. The metric is determined as the solution to an initial value problem starting at i? = and 
evolving up to -R = — rh- We indicate the location of the Gauss-Lobatto gridpoints on the spectral grid 
which is used in the pseudo-spectral collocation point method. 



which we collect in the 2n^ -dimensional vector 



Unn 



N = n-\. 



(15) 



From any such vector f the spectral coefficients and c'^ of U and Q. can be computed. Moreover, 
the corresponding coefficients of first and second derivatives of \J and i? are easily determined. 
Returning from coefficient to physical space, we can build vectors f r, . . . , f containing the values 
of first and second derivatives C/ /j, i7 /j, . . . , f/^ee, ^fiB at all collocation points 0^). 

Now consider the system (|4l|5]l at Qk) and insert the relevant components of f , f /j, . . . , f eg into 
these equations. We obtain a non-linear algebraic system of equations of the form 



F(f) = 0. 



(16) 



where 



-^00 



(17) 
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Figure 2. Initial data for Kerr and non-Kerr solution (upper panels) and corresponding values at the 
future Cauchy horizon T-L^ (lower panels). 



Here, Fj^^ and Fjj^^ denote the difference of left and right hand sides of equations ^ and ^ 
respectively, taken at the collocation point {Rj,6k). The solution f of the discrete algebraic system 
([T6l) describes the spectral approximation of the solution to our hyperbolic initial value problem. 
We find the vector f using a Newton-Raphson scheme. 



f = lim f„ , f^+i = im- [J(fm)]"' F (f^) , (18) 
where the Jacobian matrix is given by 

dF 

•>=«■ 

Note that for the convergence of the scheme a 'good' initial guess fp is necessary. For the first 
calculation, we take the corresponding coefficients of the Kerr solution for some specific parameters (say 
J and A+). We may then depart from the Kerr solution and gradually approach some new solution with 
a non-Kerr initial data set {n+(0), a;^(0)}. 

In the following example we start with a Kerr black hole in which we prescribe the specific angular 
momentum J/M^ = 0.8, where M denotes the black hole mass. As explained above, we gradually 
depart from this solution and approach an initial data set in which we replace the initial function u'^{6) 
by a constant (we take u^{6) = 0.65) but maintain uj'^ and oj~^{9) as previously determined for the 
Kerr case of J/M"^ = 0.8. The initial data u^{9) and uj^{9), as well as the functions u and u, at 
T-L~, are displayed in figure |2] Here, the values at the inner Cauchy horizon R = — rh emerge from the 
pseudo-spectral calculation. 
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Figure 3. Numerical test of relation dl} for the non-Kerr example displayed in figure 



The numerical test of relation ([T} is exhibited in figure |3] It shows a convergence plot of the expression 
|l - (87rJ)2/(A+A-)| for our sample non-Kerr solution. We obtain a confirmation of the validity of ^ 
by 12 relevant decimal digits which is equivalent to the numerical round-off error 

3. Analytical studies 

This section about the analytical study of the interior hyperbolic black hole region, including the rigorous 
derivation of the equality ([T|l, is a brief summary of earlier work that has been presented in [IJ. 

3.1. Weyl coordinates 

As a first step towards a strict mathematical treatment in terms of Bdcklund transformations, we introduce 
canonical Weyl coordinates (p, 99, t) in a small exterior vacuum vicinity of the black hole|l: 

p2 = 4(^2 _ ^2^ g.j^2 C = 2Rcose. (20) 

The corresponding line element reads as follows: 



e2'=(dp2 + dC2) + p2d(^2 _ (dt + advp)^ (21) 



Here, the metric potentials U, k and a are functions of p and The rotation axis (sin = 0) is given by: 

p = 0, |Cl>2rh, (22) 
and the exterior event horizon Ti^ is located at 

p = 0, -2rh < C < 2rh . (23) 

The event horizon is a degenerate surface in Weyl coordinates. As the interior region is characterized 
by < rh, the corresponding Weyl coordinate p assumes imaginary values there. This means that the 
hyperbolic black hole region cannot be accessed in Weyl coordinates (see figure |4l). 



^ We assume that for physically reasonable types of matter surrounding our stationary black hole, the immediate vicinity of 
the event horizon must be vacuum, see. e. g. discussion in |31 
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Figure 4. Portion of a black hole space-time in Weyl coordinates (left panel) and Boyei-Lindquist-type 
coordinates (right panel) (figure taken from |IT]). 



3.2. The Ernst equation 

The complex Ernst potential / combines the metric functions U and b, 

f = e^^ + ib, (24) 
where the twist potential b is related to the coefficient a via the relations: 

ap = pa'^^^b^c: o,C = -pe"'^^6,p. (25) 
In Boyer-Lindquist-type coordinates, this relation reads as: 

a^R = -2 sin 6 e'^^b^, a,e = 2{R^ - rg) sin 9 e'^^b^R. (26) 

The vacuum Einstein equations (which are valid in a vicinity of "H^) are equivalent to the Ernst equation 
in which reads in Weyl coordinates as 

m) {f,PP + /,CC + = /,p + fl (27) 

and in Boyer-Lindquist-type coordinates: 

(3?/) - rDf^RR + 2Rf,R + f^ee + cot Of^] = {R^ - rl)f% + f% (28) 

Note that because of the degeneracy of T-L^ in Weyl coordinates, the potential / is, for p = 0, only a 
C'^-function in terms of Q. In contrast, for a regular black hole, / is analytic with respect to the Boyer- 
Lindquist- type-coordinates R and 9. 

3.3. Bdcklund transformation 

The Ernst equation is, when written in Weyl coordinates, the integrability condition of an associated 
linear matrix problem \^ This is the great advantage of the Weyl coordinates with respect to the 
Boyer-Lindquist-type coordinates. The existence of this linear problem enables us to apply methods 
known from soliton theory. In this contribution, we are particularly interested in the so-called Backlund 
transformations through which new solutions from previously known ones are created [8 |-[12 |. As an 
example, the Kerr solution describing a rotating black hole in vacuum can be constructed from the fiat 
Minkowski metric in this manner (see e. g. lITOl ). 

In the following we use the Backlund transformation technique in order to write an arbitrary regular 
axisymmetric, stationary black hole solution / in terms of an auxiliary 'seed' potential /q. Here, /o 
describes a space-time without a black hole but with a completely regular central vacuum region. More 
specifically, /o is characterized by the following properties: 
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(i) /o is defined in a vicinity of tfie axis section p = 0, |C| < 2rh. 

(ii) In tliis vicinity, /o is an analytic function of p and C, and an even function of p. 

(iii) For /j = and |C| < 2rh, the values of /q in terms of the event horizon values of / are given by 

i [2rh(6+ + 6+) - (6+ - 6+)C] / + 4rh6+6+ 



/o 



4rh/-i[2rh(6+ + 6+) + (6+-6+)C] 



(29) 



where = 6(p = 0, C = 2rh) and 6g = h{p = 0, C = — 2rh) denote the twist potential values at 
north and south poles ofH^ . 

Now, from this Ernst potential /o the original potential / can be recovered by means of an appropriate 
Backlund transformation of the following formo: 



/o 
/o 
/o 



1 

aiAi 
A? 



1 

02^2 
^2 



1 

aiAi 
A? 



1 

02 A2 
-^2 



(30) 



where 



A, 



iz 



j = l,2, 



Ki = -2rh, 



i^2 = 2rh 



(31) 



with the complex coordinates z = p + iC, ^ 

aj^z = -{XjOj+ai 
1 



p — iC, and qi, 02 are solutions to the Riccati equations 

, fo,z , , , > N fo,z 



+ ("i + Aj)- , 
/o + /o 



— a ■ + Uj 



/0 + /0 



+ a,- + 



1 



/o,5 



A j / /o + /o 



(32) 
(33) 

(34) 

For a regular black hole, / is analytic with respect to R and cos 9 in an exterior vicinity of T-L^ . Hence, 
we can expand it analytically into an interior vicinity of 7^+. A theorem by Chrusciel (theorem 6.3 in 
lfT3l ) assures that the potential / exists as a regular solution of the interior Ernst equation for all values 



with 



1. 



iR,cose) e (-rh,rh] x [-1,1]. 



(35) 



This region only excludes the Cauchy horizon T-L^ (R = — rh). 

In the following derivation of the expression for / at the interior boundary R = — rh, a crucial role is 
played by the fact that /o is even in p. In terms of the Boyer-Lindquist-type coordinates, this means that 



/o is an analytic function of {R 



sin 6* and RcosO. The analytic expansion of /o into the region 



R < r\i retains this property. Hence, /o taken at the boundaries of the inner hyperbolic region, can be 
expressed in terms of /o taken at ii = rh. Specifically we obtain 



fo{R= - rh , cos 0) = fo{R = + rh , 



cos I 



/o(p = 0,C 



-2rh cos I 



(36) 



Also, it follows that /o is regularly defined in a sufficiently small vicinity of the boundary of the interior 
region, see figure |5] 

From the values of /o at these boundaries we can construct f onT-L^ via the Backlund transformation, 
which is stated in the following. 

A bar denotes complex conjugation. 

We obtain Chrusciel's form of the line element by substituting R = rh cos T and 9 — ip. 
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Figure 5. The seed function /q can be regularly defined in at least the grey areas (figure taken from 1 1]). 



Theorem 3.1. Any Ernst potential f of a regular, axisymmetric, and stationary black hole space-time 
with angular momentum J ^ 0, can be regularly extended into the interior of the black hole up to 
and including an interior Cauchy horizon, described by R = — in Boyer-Lindquist-type coordinates 
{R, 6). The values of f on the Cauchy horizon are given by 

f{R = -rh, cos 9) = — — (37) 

2fo{R = rh, - cos 6*) - 1 [di + 02 + (oi - 62) cos d\ 



with 



5+ - 6+ + 2(6 ' 5+ _ 6+ + 2(6 ee)^ ^^^^ 

where the scripts '+' and 'N/S' indicate that the corresponding value ofb or its second 6 -derivative has 
to be taken at the event horizon's north or south pole respectively. The values of the seed solution fofor 
R = follow via fl29D from f on the event horizon. 

Note that, for J — 0, the Cauchy horizon becomes singular. In this case we have \6i/2\ — )■ 00 and 
consequently /|^- diverges. 

By virtue of the above theorem we are able to write the inner Cauchy horizon area completely 
in terms of Ernst expressions taken at the event horizon T-L^. As discussed in 1 1 1, the universal equality 
([T]l arises in this manner. Moreover, the vanishing of is obtained in the limit J — 0, i. e. when the 
Cauchy horizon becomes singular. 



4. Discussion 

In this contribution we discussed the interior hyperbolic region of axisymmetric and stationary black 
holes which are surrounded by a matter distribution. We first looked at the corresponding degenerate 
hyperbolic Einstein equations in terms of fully pseudo-spectral methods, and confirmed the validity of 
relation ([1} to high precision. 

In the second part of the article, we used the Backlund technique in order to write the black hole 
metric in terms of an auxiliary seed potential /o that does not describe a black hole but a completely 
regular central region. We first derived simple relations between the values of /o at the boundaries of 
the interior region, and then carried these relations over to the original black hole metric by virtue of an 
appropriate Backlund transformation. A particular consequence of this relation is the universal equality 

(ID- 
Note that the key point for this method to work is the fact that the two horizons T-L^ are connected by 
specific axis sections. Indeed, the linear matrix problem associated to the Ernst equation (see discussion 
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at the beginning of section 3.3) simplifies considerably for p = 0, that is, on the two horizons Ti^, and 
on the rotation axis (sin 9 = 0). An alternative derivation, analyzing the linear matrix problem directly 
at the several sections where p = 0, was carried out in |[T4] [151 for the Einstein-Maxwell field. In that 
work, the corresponding steps yield the more general formula 

{Sirjf + {A7rQ^f = A+A- 

where Q is the electric charge of the black hole. 

It is interesting to remark that there are associated inequalities relating angular momentum and 
area of the black hole. In LI 6.1 (see also [17]) it was shown that for subextremal black holes (which 
possess trapped surfaces in every sufficiently small interior vicinity of the event horizon), the following 
inequalities hold: 

A- < v/(87rJ)2 + (4^Q2)2 < 

In the case of pure Einsteinian gravity, this relation was proved in a different context, namely for non- 
stationary, axisymmetric, vacuum space-times fTSl. In particular, the local inequality A > 87r| J| was 
shown, where A and J are the area and angular momentum of any axially symmetric closed stable 
minimal surface in an axially symmetric maximal initial data. 
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